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Abstract 

We prove that the number of limit cycles, which bifurcate from a 
two-saddle loop of an analytic plane vector field Xq, under an arbitrary 
finite-parameter analytic deformation X\, A € (1^,0), is uniformly 
bounded with respect to A. 
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1 Introduction 



Consider a finite-parameter analytic family of analytic plane vector fields 

X x = P(x,y,\)— + Q(x,y,\)—, A G R N (1) 

such that Xq has a limit periodic set T. The cyclicity of T is, roughly speak- 
ing, the maximal number of limit cycles of X\ which tend to T as A — > 0. The 
Roussarie's finite cyclicity conjecture claims that every limit periodic set oc- 
curring in an analytic finite-parameter family of plane analytic vector fields, 
has a finite cyclicity |23j . If true, the conjecture would imply the finitness 
of the maximal number H(n) of the limit cycles, which a plane polynomial 
vector field of degree n can have. Therefore it plays a fundamental role in 
all questions related to the second part of the 16th Hilbert problem and its 
ramifications. 

Recall that a polycycle of the vector field Xo is a topological polygon 
composed of separatrices and singular points. A A;-saddle cycle of Xo (or 
a hyperbolic /c-graphic) denoted Tk, is a polycycle composed of k distinct 
saddle- type singular points pi,p 2 , ■ ■ ■ ,Pk, Pk+i = Pi and separatrices (het- 
eroclinic orbits) connecting pi to p i+ i as on fig. [I] ^-saddle cycles, period 




Figure 1: One, two and three-saddle cycles. 



orbits and weak foci or centers are the simplest limit periodic sets. The finite 
cyclicity of period orbits and weak foci is well known and follows from the 
Gabrielov theorem [23, p. 68]. The finite cyclicity of one-saddle loops is due 
to Roussarie [2T1I22] . 
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The purpose of the present paper is to prove the finite cyclicity of a two- 
saddle cycle, under finite-parameter analytic deformation, see Theorem [4} 

Several special cases of this result were earlier proved, under different 
genericity assumptionss either on X or on the family X\, by Cherkas, Mour- 
tada, El Morsalani, Dumortier, Roussarie, Rousseau, Jebrane, Zoladek,Li, 
Caubergh, Luca and others 021 El El EE], see a l so [221 section 5.4.1] for 
survey of the results and references up to 1996. The finite cyclicity of a 
k-saddle cycle (any k), under finite-parameter analytic deformation was re- 
cently announced by Mourtada |18| . 

For generic families of vector fields the analitycity can be relaxed. As 
it is shown by Ilyashenko and Yakovenko |llj . Kaloshin [12] any nontrivial 
elementary poly cycle occurring in a generic k - parameter family ofC°° vector 
fields has finite cyclicity. 

In contrast to the aforementioned papers we shall not use the asymptotic 
expansions of the corresponding Dulac maps. Instead of this, we evaluate the 
number of the limit cycles near T 2 in a complex domain, by making use of 
a suitable version of the argument principle. This approach was initiated by 
the author in [10] . where we studied cyclicity of Hamiltonian two- loops. As 
it is well known, the limit cycles of planar systems close to Hamiltonian are 
closely related to the zeros of associated Abelian integrals depending on a 
parameter (the so called weakened 16th Hilbert problem [TJ Arnold, p.313]). 
Zeros of complete elliptic integrals were successfully studied by topological 
arguments in a complex domain (the argument principle) after the pioneering 
work by G. S. Petrov [191 120], see also Zoladek [25, section 6] for a description 
of the method. It has been used in a more general context in several papers, 
e.g. [31 12]; and in [10] the idea has been used to replace Abelian integrals by 
the true Poincare return map. 

In the present paper we shall find a relation between the fixed points 
of the Poincare first return map and the fixed points of holonomies of the 
separatrices of the saddle points, which correspond to complex limit cycles. 
To count such fixed points is a question on the zeros of families of analytic 
functions which is easily solved. The main technical tool is Lemma|2]in which 
we show that the connected components of the zero locus of the imaginary 
part of a Dulac map are smooth semianalytic curves. This allows to estimate 
the variation of the argument of the displacement map along the border of 
an appropriate complex domain, and finally to apply the argument principle 
in order to evaluate its zeros in the domain. 

Note that previously the relation between the monodromy and the Dulac 
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map was used by Roussarie to compute the Bautin ideal associated to the 
Dulac map [24] . This combined with [2TI |2"2"] also implies the finite cyclicity 
of one-saddle cycles. 

The paper is organized as follows. In section [2] we provide the necessary 
technical background, and prove the main technical Lemma [2j In section |3j 
we give a new self-contained proof of Roussarie's theorem about the finite 
cyclicity of one-saddle cycles. The origin of our method is then explained 
in section [4] where we give a brief account of a local version of the so called 
"Petrov trick". The same method is then easily adapted in section [5] to show 
that the cyclicity of 1?2 is finite. 

2 The Dulac map 

Consider an analytic family of plane real analytic foliations F\, A G K , 
having a non-degenerate isolated saddle point. An appropriate translation 
analytically depending on A will place the saddle point at the origin. The 
foliation F\ has two analytic separatrices, transversally intersecting at the 
saddle point, and depending analytically on A (5J HTj. Therefore a further 
real bi-analitic change of the variables x, y, analytically depending on A, will 
identify them to the axes {x = 0} and {y = 0} as on fig|2j so 

F x : ar(l + ...)dy + a(\)y(l + ...)dx, a(0) > 0. (2) 

where the dots replace higher order terms in x, y with coefficients depending 
on A. The number a(A) is the hyperbolic ratio of the saddle point. From now 
on we shall suppose that the foliation (|2j is analytic and depends analytically 
in A in a neighborhood of the origin in M 2 x M. N . 

For C\,C2 G K sufficiently small, let a C {?/ = cj, r C {x = C2} be 
open complex discs centered at (0, C\) and (C2, 0), parameterized by x and y 
respectively. The (real) Dulac map is the germ of analytic map at x = 

P A :ffnR + ^rnl + , £> A (0)=0 

defined as follows: if x G crflM^ then T>\(x) G rflMj" is the intersection with 
rfllR^ of the orbit r y\(x) of ([2]), passing through x, see|2](ii). This geometric 
definition of T>\ allows to control to a certain extent its analytic continuation 
in a complex domain. 
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(i) (H) 

Figure 2: The Dulac map 
2.1 Analytic continuation 

The Dulac map allows an analytic continuation on some open subset of the 
universal covering space a, of a \ {0}, depending on A. Let us parameterize 
cr. by polar coordinates p > 0, <p G R, z = pexpiip. The following result is 
well known (e.g. |10l Appendix A]) 

Theorem 1 There exists Eq > and a continuous function 

suc/i i/iai i/ie Dulac map allows an analytic continuation in the domain 

{(A, p, </?) G x a. : |A| < £ , < p < p(cp)} (3) 

The geometric content of Theorem [T] is as follows. Let {'jx(z)} Zt x be a 
continuous family of paths contained in the leaves of F\, and connecting 
2 G cr to T. 

For z G a R we suppose that 7^ (z) is the real orbit of JF\ contained 
in the first quadrant x > 0,y > 0, and connecting z to r, see fig{2] (ii). 
The above Theorem claims that this family of orbits allows an extension to 
a continuous family of paths {7a(z)},z,a, contained in the leaves of F\, and 
connecting z G a. to r.. The family is defined for all (A,p, <p) which belong 
to the domain (|3]). Each path starts at z and terminates at a unique point 
on a, denoted T>\(z). Although the paths {7a(-2)} 2 ,a are not unique, their 
relative homotopy classes are uniquely defined. 
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2.2 Monodromy of the Dulac map and holonomy of sep- 
aratrices 

To the axes {x = 0}, {y = 0} parameterized by y and x, we associate holon- 
omy maps 

h\ : cr ->■ a, : r -> r 

defined by two closed paths contained in the axes {x = 0} and {y = 0} 
and based at (0,Ci), (c2,0) respectively. We shall make the convention, that 
each closed path makes one turn around the origin of the axe in which it is 
contained, in a positive direction (recall that the axes are parameterized by 
y and x respectively). It is easily seen that in the case of a linear foliation of 
the form 

xdy + a ydx = 0, a G M + (4) 

we have 

V a : x i-> y = Cl Cz a x a , h a : x xe~ 2ni/a , h T : y ^ ye" 2 ™. (5) 

In the general case of a nonlinear foliation of the form ^ the Dulac map 
D\ is only asymptotic to a , while the holonomy maps are analytic in 

x, y, A and 

xe- 2 " i/Q + . . . , : y ^ ye" 2 ™ + ...,« = a(A). (6) 

The Dulac map T>\ is a transcendental multi- valued map. For x > let 
^A(e 27rl x) be the result of the analytic continuation of T>\ along an arc of 
radius x and angle 2iri. Similarly, for y > let D\(e 2m y) be the result of the 
analytic continuation of T>\ along an arc of radius y and angle 

Lemma 1 For every sufficiently small x > 0, y > 0, |A| holds 

h x T o V x (e 2m x) = V x (x), h x a o Vl\e 2m y) = V^(y). 

Proof. Consider, instead of T>\ the underlying path 7^. The loop r )\{e 2 ^ l x) 
has the same origin as ^\{x) so they can be composed and the resulting loop 
7a (y) starts at y = V\(e 2m x) G r and terminates at T>\(x) G r. In the special 
linear case ^ with a = 1 the foliation is a fibration, the paths 7a(-) represent 
relative cycles in the fibers oixy, and the path j\(y) is closed and represents a 
vanishing cycle. The claim of Lemma [T] is then the classical Picard-Lefschetz 
formula. In the general case the result follows "by deformation". Indeed, 
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in the linear case Q with a = 1 the family of closed paths {7^ (y)} y is 
defined for all sufficiently small y, and 7a (0) C {y = 0} is a closed path 
which makes one turn around the origin in the axe {y = 0} in a positive 
direction. Note that the paths 7a (0) are bounded away from the origin in 
C 2 . It follows that 7^ defines the holonomy h x of the separatrix {y = 0}, 
and this property holds true also for every sufficiently small deformation of 
(|4]). The homothety (x,y) — > (ex,ey) transforms ^ to a small deformation 
of ^ which completes the proof of the first identity (but see also [IS]). The 
second identity in Lemma [T] is proved in a similar way. □ 

2.3 The zero locus of the imaginary part of the Dulac 
map 

Consider the universal covering 

C,4C \ {0} (7) 

and the zero locus %x C C. of the imaginary part of the Dulac map T>\ 
corresponding to the domain ^ 

H x = {z = (p, p) G C. : lmV x (z) = 0, < p < p(<p), p G R}. (8) 

In the case of a linear foliation Q the zero locus is therefore a union of 
half-lines: 

A/7T 

Ha = {z G C. : Im z a = 0} = Uftezfta,*, U a>k = {(p, <p) G C. : p = — }. 

a 

To describe H\ in the case of a general foliation of the form (|2]), with hyper- 
bolic ratio a (A) > 0, consider the germs of real analytic sets at the origin in 

R 2 = C 

C x , k = {zeC = R 2 : (h x a ) k (z) = z] (9) 
where h x is the holonomy map associated to the separatrix {x = 0}. 

Lemma 2 The zero locus Tix C C, of the imaginary part of the Dulac map 
in the domain (f51) is a union of connected components Hx^, indexed by k G Z. 

• Each set Cx,k, is a germ of a real analytic curve of R 2 , which is 
smooth at the origin and tangent to the line 

{z = se ik7T/a{x) : s E (R,0)} (10) 

there. 
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Each connected component 7ix,k ^ s projected on the plane C = IR 2 under 
the map ir _{7p to the connected component of \ {0} tangent to the 
half-line (10), s > 0, at the origin. 



Remark 1 For a general bi-holomorphic map h^, vanishing at the origin, the 
set ([9]) coincides with the origin itself. The above Lemma shows, however, 
that for the monodromy map of a saddle point of a real- analytic plane 
vector field, the set C\^, is a germ of a real analytic curve ofM 2 , which 
is smooth at the origin. The position of the connected components ofC\^\{0} 
tangent to the half-lines (10), s > (10), s > is shown on fig. [31 




Figure 3: The zero locus Ti\ of the imaginary part of the Dulac map, pro- 
jected on the complex plane C. 



The above Lemma is the main technical result of the present paper. The 
analyticity of the zero locus 7-L\ is responsible for the algebraic-like behavior 
of the Dulac map. 

Proof of Lemma @ Let x G a R 1R + and suppose that for some ip > 0, 
T>\(e l(p x) E M. As the Dulac map is real along a PI M + , then V\{e~ lip x) is 
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complex conjugate to T>\(e llf x) and hence 



V x (e-^x)=V x (e^x). 

If the point e~ lip x is seen as the inverse image of T>\(e~ lip x) with respect to the 
Dulac map T>^ , then the point e ltp x is the result of the analytic continuation 
of the map T>^ 1 along a suitable closed path of r, starting and terminating 
at V x (e~ iip x). If we put 

y = V x (e-^x), e- l *x = V-\y) 

then e ±lip x are two values of the multivalued map V^iy) and hence, by 
Lemma [TJ they differ by a power of the monodromy 

(h£) k {e itp x) = e~^x 

or equivalently 

(h^) k (z) — z, z = e iip x, for some k <E Z. 

Clearly every such relation will correspond to a connected component Hx^ of 
%\. As H\ t k is an analytic set of real dimension one, then is an analytic 
set of dimension one too. It can be defined therefore by each of the following 
equivalent relations 

C x>k G{zeC = R 2 : Re[(h x a ) k (z)} = Re(z)} 



or 



C x>k C {z G C = I lm[(h x ) k (z)] = Im(z)}. 



As ■§=[(h x .) k (z) — z] = —1, then the linear part of the complex-analytic func- 
tion 

R 2 -»■ C 
(z,z) ^ (ti) k (z)-z 

can not be identically zero, and therefore C\j~ C 1R 2 is a real analytic curve, 
smooth at the origin. It follows from ^ that the projection of Hx,k under 7r 



on the plane C = R is tangent to the half-line (10), s > at the origin. □ 
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2.4 The argument principle 

Let D C C be a relatively compact domain, with piece-wise smooth bound- 
ary, and ip : D — > C an analytic function which allows a continuation to the 
closure D. Denote by Z-£>(ijj) the number of the zeros of ip in D, counted 
with multiplicity. If we assume that ip does not vanish along the border <9D, 
then the increment of the argument VardD{aig(ip)) of ij) along <9D oriented 
counter-clockwise is well defined. V ar QY)(axg(ip)) equals the winding number 
of the curve ijj(dD) C C about the origin and the classical argument principle 
states that 

2vrZ D (^) = ^ar 9D (arg(^)). (11) 

More generally, if ip has zeros on <9D, isolated or not, the variation of the 
argument V ar QY>(sxg(%p)) might be not well defined. 

Definition 1 We say that z G <9D is a regular zero of ip if tp(z) = 0, and i[> 
allows an analytic continuation in a neighborhood of z in C 

If we assume that ip has only regular zeros in D, then VarQn(axg(ifj)) is 
well defined as a sum of the increments of the argument of ip\dr> between 
consecutive zeros of ip. Indeed, the increments are finite, because the border 
dD is piece-wise smooth. The argument principle can be reformulated as 
follows 

Proposition 1 Let D C C be a relatively compact domain with piece-wise 
smooth boundary. Ifip : D — > C is a continuous function, analytic in D ; and 
having only regular zeros in D, then 

2ttZ d (^) < Var dD (arg(^)) < 2ttZ d (^) + 2ttZ 9d (^) (12) 

Proof. There always exists a polynomial P, such that ip/P has no zeros 



in D, so we need to verify (12) for polynomials only. The set D is open, 
connected and oriented, with piece-wise smooth boundary, which therefore 
has no self-intersections and has an induced orientation. The inequality 

< Var dD (arg(z)) < 2vr 



allows to "remove" the zeros along dD and hence formula (11) implies (12) 
□ 



In the present paper the first inequality in (12) will be used to bound the 



number of the zeros Zd(.). For this we shall need estimates on the variation 
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of the argument V ari(arg(.)) along any compact segment I of a curve. More 
precisely, let I C R 2 = C be a compact segment of a smooth real analytic 
curve. Let U C C be an open set containing I and ip\(z), A G (C N ,0), be 
a germ of a family of complex-analytic functions in U at A = 0. For every 
fixed A such that the function ip\ is not identically zero, the variation of its 
argument 

|Vari(arg(V>A)| 

is well defined. 

Theorem 2 Let I be a compact segment of a real analytic curve and let 
{^a}a be a family of functions analytic in a neighborhood of I, and depending 
analytically in A. There exists Eq > 0, such that 

sup \Vari(axg(ip\)\ < oo. 
|A|<£o,Va^o 

The above result follows from the following theorem due to Gabrielov |14| E] 

Theorem 3 Let M, N be real analytic varieties and consider the canonical 
projection tt : M x N — > N. For every relatively compact semianalytic set 
E C M x N , the number of the connected components of the pre-images 
7r _1 (ra) is bounded from above uniformly over n G N. 

Proof of Theorem [2] . The number of the isolated zeros of ip\ along I 
counted with multiplicity is uniformly bounded in A at A = (Frangoise- 
Yomdin Theorem |14j). On an interval between two zeros of the varia- 
tion of the argument divided by 2ir is bounded by the number of the zeros of 
the imaginary part of ip\ divided by two, plus the sum of the multiplicities 
of the zeros of tp\ at the end of the interval. The imaginary part of ipx is 
a real analytic function in U C 1R 2 and the Gabrielov Theorem implies that 
the number of the connected components of {Im(ip\) = 0} D I is uniformly 
bounded in A at A = 0.D 

3 Cyclicity of one-saddle cycles 

Let Xx, A G (M. N , 0) be a germ of an analytic family of analytic plane vector 
fields, such that Xq has a one-saddle cycle (homoclinic saddle loop) Vi. The 
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first-return map associated to Fx is a composition of a Dulac map T>\(z) : 
a — > r and a transport map T\{z), see fig. |4j We assume that the Dulac 
map is in a normal form as in section |2.1| The limit cycles of X\ near I\ 
correspond to the zeros of the displacement map 

i} X (z)=V x (z)-Tx(z) 

near z = 0. An appropriate choice of the local coordinates on the croos- 
sections a and r brings the transport map to the form 7~\(z) = z. Alter- 
natively, we could choose simply a = r (without supposing that the Dulac 



map is in the normal form of section 2.1). We shall bound the number of 



the zeros of ip\ in the domain D^j C C. delimited by the circle {p = R}, and 
the connected components H\\ and "Ha-i of the zero locus of the imaginary 
part of the Dulac map, as it is shown on fig. [5j We shall suppose that R > 
is so small, that is analytic in D# for all A G Mr, such that |A| < Eq 

(Theorem [T]) and it is analytic even on the closure of D# except of course at 
z — 0, where ipx(-) is only continuous. Indeed, 

lim V x (z) = 

z^0,zeD]j 

while T\(z) is holomorphic at z — 0, so 

lim i[)x(z) = c(A) 
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Figure 5: Examples of domains D# C C., projected on the complex plane C 
under it 

where c(A) is analytic and c(0) = 0. If the family of functions ip\ is sufficiently 
general, then c(A) ^ 0, and in the case when c(A) = we may replace ip\ 
by the new family ip\ + Ajv+i, A^r+i G K, for which the limit at z = is the 
parameter Ajv+i- After this preparation, we may prove the finite cyclicity 
of the homoclinic loop Ti. For this we apply Proposition [T] (the argument 
principle) to the family of functions ip\ in the domain D#. In the course of 
the computation, it will be supposed that R > is sufficiently small, Eq is 
sufficiently small with respect to R, and A is such that |A| < Sq. We may 
encode this choice of the parameters by the "physical" notation 

< |A| < £ « R « 1. (13) 

The hyperbolic ratio of the saddle point will be not bigger than one only 
in a suitable semi-analytic set in the parameter space, and will be bigger than 
one in another (complementary) semi-analytic set. After eventual exchanging 
of a and r, it will be also supposed that the hyperbolic ratio of the saddle 
point is not bigger than one for all parameter values. 

Along the circle {z : \z\ = R} with angle close or strictly less than 2ir the 
variation of the argument of ip\ is uniformly bounded in A (Theorem pi). 

Along the curve C\,i the imaginary part of ip\ equals the imaginary part 
of the transport map — T\(z) = —z. Therefore the zeros of Im(tfj\) along 
C*a,i are exactly the intersection points of C\,i and the segment (— R,0). 
According to Lemma [2] we have 

C x ,i DR = {x ER: h*(x) = x} = C x ,-i D K. (14) 
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As h*(x) is an analytic family of analytic functions, then by Gabrielov's the- 
orem, the number of such fixed points is uniformly bounded in A on [-R, 0] . 
To conclude, we have only to check that the family has regular zeros 

along the border of the domain D#. This is indeed the case, when c(A) 7^ 0, 
as ip\(0) = c(A). We conclude that the number of isolated zeros of the family 
of functions 

{^ A :c(A)^0,|A|<£ } 

in the domain is uniformly bounded by some integer, say C. Finally, 
note that the condition c(A) 7^ can be removed. Indeed, if for some Ao, 
|Ao| < ^O) c (-^) = 0, the function ip\ has at least C + 1 zeros in D^j, then it 
has at least C + 1 zeros in Dr in a sufficiently small neighborhood of Ao, in 
contradiction with the preceding estimate. 

To resume, we proved the following classical result 
Theorem (Roussarie|21], 1221 |24| ) Every homoclinic saddle loop (a one- 
saddle cycle) occurring in an analytic finite-parameter family of plane ana- 
lytic vector fields, may generate no more than a finite number of limit cycles 
within the family. 

Let us note that our method, exactly as the Roussarie's Theorem allows 
to compute more precisely the cyclicity of I\. We shall not enter into details 
here. Just to illustrate this, note that if the hyperbolic ratio ct(0) is strictly 
bigger than one, then the overall increase of the argument of the displacement 
map along the border of is strictly less than 2it (this computation is 
omitted) and the cyclicity of Ti is zero. 

4 The Petrov trick 

The content of this section is not necessary for the proof of our main result 
Theorem [4], but it aims to shed some light on the origin of the method, used 
to bound the limit cycles near the saddle loop in the preceding section. 

With the same notations as in section [3j consider the analytic family of 
analytic vector fields 

X x , X = (X 1 ,...,X N )e(R N ,0) 
defining a holomorphic foliation T\ of the form 

F x = {dH + \ 1 u x = 0},u ^0 
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where H is a function and u\ is an analytic family of differential one-forms, 
both analytic in a neighborhood of the saddle loop Y\. For definiteness, we 
put the saddle point at the origin in IR 2 , so dH(0) = 0. We shall further 
suppose that the saddle loop Ti is contained in the level set {H(x,y) = 0}, 
and the interior of r\ is filled up by a continuous family of periodic orbits 
Jo(h) C {H(x,y) = h}, parameterized by h > 0, where h = H(x,y)\ a 
is the restriction of H on the cross-section a. The displacement map is 
approximated by the usual Poincare-Pontryagin formula as follows 



where o(Ai)/Ai tends to zero as A tends to zero, uniformly in h in every 
compact interval in which the displacement map is defined. The zeros of 
ipx(-) correspond to limit cycles and, at least far from Ti C {H(x,y) = 0}, 
they are approximated there by the zeros of the complete Abelian integral 



We make the assumption (actually justified by the Roussarie's theorem |21|). 
that this is so also in a neighborhood of h — (corresponding to limit cycles 
close to the saddle loop Ti). Thus, it makes a sense to prove the finiteness of 
the maximal number of the zeros of the Abelian integral I\(h), which tend to 
h = as A tends to the origin in the parameter space. This follows of course 
from a well known general result of Varchenko and Khovansky. We shall 
use, however, a different idea due to G.S. Petrov (20] , who showed that the 
analogous global problem for complete elliptic integrals of second kind is of 
algebraic nature. This observation has been used in several papers by Petrov 
to evaluate the precise number of zeros of complete elliptic integrals, and 
hence of limit cycles of perturbed Hamiltonian vector fields, see for instance 
Zoladek section 6]. We are ready to describe the local version of the 
Petrov method. 




(15) 




Consider the sector 



S R = { Z = pe lv eC:0 <p<R, 0<y?< 2vr}. 



For a fixed sufficiently small R > and all sufficiently small ||A|| 
the Abelian integral I\ (z) allows an analytic continuation in Sr. 
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To bound the number of its zeros on Sr (and hence on (0, R)) we 
apply the argument principle to the domain Sr. Along the circle 
{p = R} the increase of the argument of I\ is bounded uniformly 
in A (due to Gabrielov's theorem). Along the segment [— R, 0] 
the Abelian integral allows two analytic continuations I^ih)- As 
I\(.) is real-analytic on (0,R) then 



I+(h) = l;(h), he(-R,0) 

and by the Picard-Lefschetz formula 
2 v / ^Im/+(/i) = /+(/>)- /-(/>)= / u x ,hE(-R,0) (16) 

J&(h) 

where 5(h) C {H(x,y) = h} is a continuous family of cycles, 
vanishing at the origin as h tends to zero. 

The imaginary part of I\(h) on (— R, 0) is therefore an an- 
alytic function, and by Gabrielov's theorem again, its zeros are 
uniformly bounded in A on the closed interval [-R, 0] . This im- 
plies that the increase of the argument of I\(h) on (— R, 0) is 
also uniformly bounded in A which combined to the argument 
principle shows the fmiteness of the maximal number of zeros 

The proof of the finite cyclicity of the one-saddle loop from the preceding 
section, may be seen as a generalization of the Petrov method. Indeed, 
the Picard-Lefschetz formula corresponds to the claim of Lemma 1, and by 



Lemma 2 the zeros the analytic Abelian integral (16) correspond to the fixed 
points (complex limit cycles) of the holonomy map h\ of the separatrix. As 
it is well known, the holonomy map of a separatrix is analytic, which implies 
the finite cyclicity of the saddle loop Y\. 



5 Cyclicity of two-saddle cycles 

The main result of the paper is the following 

Theorem 4 Every heteroclinic saddle loop (a two-saddle cycle) occurring 
in an analytic finite-parameter family of plane analytic vector fields, may 
generate no more than a finite number of limit cycles within the family. 
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Figure 6: The Dulac maps T>\ and V\ 



Using the notations used in the preceding sections, suppose that the vec- 
tor field Xo has a two-saddle loop r 2 . Consider the Dulac maps 

V\:a^T, i = l,2 

associated to the corresponding foliation, as on fig 6j Each map T>\ is a 
composition of a "local" Dulac map (as in section[2) and two real-analytic 
transport maps. From this it follows that Lemma ^applies to T>\, i = 1,2, 
too. From now on we choose a real-analytic local variable z on the cross- 
section a thus identifying a to an open disc centered at G C. We shall also 
suppose that = cr nr 2 . The functions V\(z), i = 1, 2 are multivalued on the 
cross-section a and have critical points at Sj(A) G 1R, Sj(0) = 0, respectively. 
The functions Si are real-analytic. The limit cycles of X\ near T 2 correspond 
to the zeros of the displacement map 

^ x {z) = V\{z)-Vl{z) 

near z = 0. Let «j(A) > 0, i = 1,2 be the hyperbolic ratios of the saddles. 
We shall suppose, upon exchanging eventually the roles of a and r, that 
cvi(0)a!2(0) > 1. Denote the zero loci of the imaginary parts of the Dulac 
maps D\(z), T>\(z) by H\ and respectively. We shall bound the number 
of the zeros of ip\ in the complex domain of the universal covering of 
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C \ {si (A), 32(A)} defined as follows (without loss of generality we assume 
that si(A) < s 2 (A)). 



if 0(2(0) > I? the domain D^j is bounded by the circle 

S R = {z : |z| = R}, 

and by 



(17) 



rl -A,l> A, — 1 > n \,-l 



as it is shown on figjTj 



• if 02(0) < 1 then necessarily Oi(0) > 1. The domain D^j is bounded 
by the circle Sr, by the interval [si(A), s 2 (A)], and by H\ 1 ,H\ _ x , as it 
is shown on figjSj 




n 2 



^A,-l 



Figure 7: The domain D# C C. projected on the complex plane C in the 

case 02(0) > 1. 



In the course of the proof the parameters R and A will be chosen as in 
the one-saddle case: the constant R will be sufficiently small, e$ > will be 
sufficiently small with respect to R, and A G M. N will be such that |A| < Eq, 
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see (13). Like in section |3| we shall suppose, without loss of generality, that 



the analytic functions c\ (A), 02(A) where 

lim i/; x (z) = ci(A), lim f x {z) = c 2 (A) 

z— >si(A),zSD fl z— >S2(A),zeD B 

are not identically zero. This will guarantee that for generic values of A the 
displacement map will have only regular zeros in the closure of D#, so the 
argument principle (Proposition [T]) can be applied. 

Proof of Theorem [H It follows from the definition of the domain 
Db C C, that the displacement map ip\(z) is analytic there. To count the 
zeros (corresponding to real and complex limit cycles) of the displacement 
map in D^j we apply Proposition 1 (the argument principle) to the family of 
functions ip x . To evaluate the variation of the argument of the displacement 
map along the border of D/j we repeat the arguments of section [3j 

Consider first the case 0:2(0) > 1, fig {7} The connected component of 
the zero locus of the imaginary part of T>\ which is tangent to the line ip = 
7r/o 2 (A) through s 2 (A) intersects the circle Sr transversally, and along this 
circle the variation of the argument of ip x is uniformly bounded in A (Theorem 
[2]). The imaginary part of ip x (z) restricted to 1-L\ equals the imaginary part of 
— D\ and hence Im^ vanishes along 'H\ 1 ,'H\_ 1 exactly at the intersection 
points 

According to Lemma [2] these intersection points are the solutions of the 
equation 

h&z) = h*(z) (18) 

where h\, h\ are the holonomies of the separatrices intersecting a and related 
to the saddle points Si(A) and 52(A). By Gabrielov's theorem, the number 
of such fixed points is uniformly bounded in the disc {z : \z\ < R}. 

Consider now the second case 02(0) < 1, «i(0) > 1, see figj8} Along 
this circle Sr the variation of the argument of ip x is uniformly bounded in 
A (Theorem pi). Along the interval [si(A), s 2 (A)] the imaginary part of D\ 
vanishes identically, and the imaginary part of ^(z) restricted to this interval 
equals the imaginary part of — V>\. Therefore the zeros of lm(ip x ) along 
[si(A), s 2 (A)] are exactly the intersection points of H X1 and [si(A), s 2 (A)]. 
By Lemma [2j and like in (14), these intersection points are the solution of 



the equation 

hn(z) = z 
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-1 



Figure 8: The domain D# C C, projected on the complex plane C in the 

case 0:2(0) < 1, «i(0) > 1. 



where h" is the holonomy of the separatrix intersecting a and related to the 
saddle points 52(A). By Gabrielov's theorem, the number of such fixed points 
is uniformly bounded. Finally, the zeros of Im(^^) along T-L\ 1 and T-L\ _ 1 are 
evaluated as in the case 02(0) > 1. This completes the proof of Theorem 131 □ 



6 Concluding remarks. 



The identity (18) which determines complex limit cycles "responsible" for 



the cyclicity of the double loop T 2 is the main new ingredient in the proof 



with respect to the one-saddle case. Indeed, solutions of (18) are fixed points 
of the holonomy h\ o ih\)~ x which, for A = 0, is generated by a closed loop 7 
contained in the complexified separatrix of T 2 intersecting the cross-section a. 
The topological type of this separatrix near T 2 is a disc with two punctures, 
corresponding to the two saddle points Si (A) and S^A). Clearly 7 makes 
one turn around each of them, but depending on the orientation we have two 
possibilities shown on fig. [9] (i) and (ii). A simple computation on a model 
example shows that the loop 7 associated to the holonomy ° (hi) -1 is the 
figure eight-loop on fig. [9] (i). The reader will recognize in the loop 7 a key 
ingredient in the proof of the local boundedness of the number of zeros of 
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pseudo-Abelian integrals in [3J H] . 




Although the result of Theorem [4] is existential, the proof we use leads to 
effective upper bounds on the number of the bifurcating limit cycles. This 
possibility is explored in [10], where we show that the cyclicity of a Hamilto- 
nian two-loop is bounded by the number of the zeros of a pair of associated 
Abelian integrals, a phenomenon which also explains the appearance of alien 
limit cycles in |S]. 



Figure 10: Hyperbolic planar polycycles with finite cyclicity. 



It worth noting, that our finitness result holds true, with the same proof, 
for other hyperbolic polycycles (on the plane or on an analytic surface), as 
those shown on fig. [10 
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